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ABSTRACT 

The chemistry of astronomical hydrocarbons, responsible for the well-known infrared emission 
features detected in a wide variety of targets, remains enigmatic. Here we focus on the group of 
young intermediate-mass Herbig Ae stars. We have analyzed the aliphatic and polycyclic aromatic 
hydrocarbon (PAH) emission features in the infrared spectra of a sample of 53 Herbig Ae stars, 
obtained with the Infrared Spectrograph aboard the Spitzer Space Telescope. We confirm that 
the PAH-to-stellar luminosity ratio is higher in targets with a flared dust disk. However, a 
few sources with a flattened dust disk still show relatively strong PAH emission. Since PAH 
molecules trace the gas disk, this indicates that gas disks may still be flared, while the dust disk 
has settled due to grain growth. There are indications that the strength of the 11.3-^m feature 
also depends on dust disk structure, with flattened disks being less bright in this feature. We 
confirm that the CC bond features at 6.2 and 7.8 fim shift to redder wavelengths with decreasing 
stellar effective temperature. Moreover, we show that this redshift is accompanied by a relative 
increase of aliphatic CH emission and a decrease of the aromatic 8.6-/im CH feature strength. 
Cool stars in our sample are surrounded by hydrocarbons with a high aliphatic/aromatic CH 
ratio and a low aromatic CH/CC ratio, and vice versa for the hot stars. We conclude that, while 
the overall hydrocarbon emission strength depends on the dust disk's geometry, the relative 
differences seen in the IR emission features in disks around Herbig Ae stars are mainly due to 
chemical differences of the hydrocarbon molecules induced by the stellar UV field. Strong UV flux 
reduces the aliphatic component and emphasizes the spectral signature of the aromatic molecules 
in the IR spectra. 

Subject headings: astrochemistry, stars: pre-main sequence, planetary systems: protoplanetary disks, 
infrared: general 

1. Introduction 
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Herbig Ae stars are pre- and zero-age main- 
sequence stars of a few solar masses. They 
are surrounded by a circumstellar disk, a rem- 



nant o f the star-formation process. iMeeus et al.l 
( 200l[ ) have classified the spectral energy distribu- 
tions (SEDs) of these targets into two groups, 
which reflect the geometry of the dust disk. 
Group I sources have a strongly flared outer disk, 
while the disks in group II sources are flatter. 
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The outer disk surface of the latter hes in the 



shadow of the inner disk toullemond et al. 


2001; 


Duhemond & Dominik 2004a; Acke et al.l 


2om. 


Both theory (DuUemond & DominikI 2004U 


Meiier et al. 


20081) and observations ( Acke et alj 2004 ) indicate 



that this difference in disk geometry is due to dust 
grain growth and subsequent settling to the mid- 
plane. Also the smallest dust grains in the disk 
surface layer appear to show evidence for sedimen- 
tation, although the correlation between grain size 
and flaring co uld also be due t o a radi al gradient 
in gr ain size ( Bouwman et al.l 2008: F urlan et al 



20091 Juhasz et al. 2010) 



The material content of these circumstellar 
disks has been primarily studied based on in- 
frared (IR) spectra obtained with ground-based 
instruments and spectrographs aboard satellite 
missions su ch as the In f rared Space Observa- 
tory (ISO, iKessler et al. 1996[ ) ai id, more re- 



cently , the Spitzer Space Telescope ([Werner et al 
2004al ). IR spectroscopy probes the thermal 



emission of the warm surface layer of the disk. 
Micrometer-sized silicate dust grains, both amor- 
phous and crystalline, are commonly detected 
both because of their abundance and the fact 
that their optical properties give rise to strong 
ba nds in this wavelength r egime (see the review 
by iHenning fc Meeud 2009h . The spectral finger- 
print of the silicates has been interpreted in terms 
of grain growth and cryst a.llisation of the small 



est grains in the disk (e . g., iBquwman et al.ll2001 



Acke fc van den Ancker 



20041: 



van Boekel et al 



2005t ISargent et al. 2009). Besides silicate bands. 



the IR spectra of many Herbig Ae stars display 
strong emission features from polycyclic aromatic 
hydrocarbon (PAH) molecules. In this paper, we 
focus on these molecules. They trace the carbona- 
ceous component of the circumstellar matter. 

Next to its detection in Herbig Ae stars, 
PAH emission is observed in a multitude of ob- 
jects and environments, including the interstellar 
medium (ISM), T Tauri stars, reflection nebu- 
lae (RNe), and, occasionally, evolved stars. The 
emission features are the result of IR fluorescence 
of PAH molecules pumped by ultraviolet radia- 
tion. Ionized PAHs can also be excited by opti- 
cal and near-IR photo ns (e.g.. iLi fc Drainel 2002 : 
Mattioda et al. 2005al lb). Because of the stochas- 
tic nature of the heating process, PAH emission 
can be seen even far from the illuminating UV 



source. For a compr ehensive review on interstel- 
lar PAHs, we refer to lTielensI (|2008l ). 

Astronomical PAH features come in different 
flavors and show a broad variety in shape, wave- 
length position and relative intensity. It is be- 
lieved that this diversity has its origin in a dif- 
ferent chemical constitution and ionization state 
of the carriers. The wealth of laboratory mea- 
surements and theoretical computations of hydro- 
carbon IR spectra has deepened our understand- 
ing of the emission features. While the features 
at 3.3 /im, 8.6 /im and in the 10-14-/im range are 
attributed to CH stretching and bending modes, 
the features at 6.2 /im a nd 7-9 Atm are due to CC 
stretching modes (lAllamandola et al. 19891) . Ion- 
ization enhances the strength of the CC modes 
(6-9 /im) relative to the CH modes (10-14 /im). 
The molecular size mainly affects the CC fea- 
tures, while the changes in the CH modes are ex- 
pected to be more subtle. The 8.6-/^m feature, at- 
tributed to the CH in-plane bending mode, is sug- 
gested to be prod uced by large (N,^^ 100 ) com- 
pact PAHs only teausc hlicher et al.llioo i). The 
7-9- fim complex consists of several sub-bands, 
with the most prominent at 7.6 and 7.8 /im pro- 
duced by small and large ionized PAHs respec- 
tively (|Bauschlicher et al. 2008, 2009). 

Despite our theoretical and experimental knowl- 
edge of PAHs, none of the astronomical spectra 
can be satisfactorily fitted with the available set 
of computed/laboratory spectra. Many analyses 
therefore focus on observational correlations. Fea- 
tures produced by CC modes correlate better with 
each ot her than with CH- mode features, and vice 
versa. Honv et "aD ( 200l[ ) have investigated the 



CH out-of-plane modes (10-14 /j,m) in a sample 
of various astronomical targets. The 11.3-//m fea- 
ture is ascribed to lone (solo) CH groups found in 
PAH molecules with long straight edges, while the 
features at 12.0, 12.7, 13.5 and 14.2 ^m are linked 
to duo, trio, quartet and quintet CH modes, in- 
dicative of corners in the molecular structure. 



Peeters et al.l (|2002l ) categorized the astronom- 
ical PAH spectra in 3 classes, reflecting the peak 
position of the major features in the 6-9- /im range. 
Class A contains the targets with the bluest fea- 
tures, mainly general ISM sources, Hii regions 
and RNe. Class C, characterized by the red- 
dest features, consists mainly o f cool carbon-rich 
(post-)AGB and red giant stars ( Sloan et al.ll2007 : 
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Gielen et all 120091: ISmolders et"alll2010l ): Class B 
sources have intermediate feature peak positions, 
and include theHerbig Ae st a rs and most PNe. 
Sloan et alj (|2005[) and ISloanI (|2008[) propose to 



categorize the PAH spectra of Herbig Ae stars in 
a new class (dubbed Class B'). Their 7-9-^m fea- 
tures show an extended red wing beyond 8.0 /Ltm, 
which is absent in the Class B features of PNe. 

The wavelength shift of the features has been 
observationally lin ked to the effective temperature 
of the central star (Slo an et al .120071: iBoersma et al. 
[2OO8; .Kcllcr_ct al 200i), with cooler stars display- 
ing redshifted IR emission features. This indicates 
that the chemistry of the emitting hydrocarbon 
molecules depends on the properties of the stellar 
radiation field. 

A plateau with superimposed narrow fea- 
tures has been detected i n the 15-21-Mni range 



in a number of targ e ts (IBeintema et al. 



Van Kerckhovcn et al. '2000*; 'Sturm et al 



1996; 



200c ; 



Werner et alll2004b:,SeUgren et al..,2007i) . It is our 



observation that m ost of these targets, if not all, 
belong to Class A. iPeeters et al.l (|2004l ) attribute 
the narrow features to out-of-plane skeletal CCC 
modes of gas-phase PAH molecules. The authors 
argue that the latter are large elongated PAHs 
with straight edges and few pendant rings. 

All features mentioned above are ascribed 
to PAHs. Their fundamental constituents are 
aromatic (i.e. benzenoid) rings. Aliphatic hy- 
drocarbon molecules do not have such rings. 
The features at 6.8 and 7.2 /xm, detected in a 
number of sources, are attrib uted to aliphatic 

Chiar et al.ll2000l: 



19991 : 



bonds (e.g.. iFurton et al 

Dartois et all 20071 : however, see also Bauschlicher et 



of shadowing of the outer disk and a change in the 
shape of the infrared continuum emission. Juhasz 
et al. (2010, submitted to ApJ; hereafter JIO) 
present the full spectral decomposition into dust 
and continuum components. The characteristics 
of the amorphous and crystalline silicate content 
of the upper disk layers are related to stellar and 
disk properties. The current paper focuses on 
the IR emission features produced by aromatic 
and aliphatic hydrocarbon molecules. Given their 
importance in the surface layers of circumstellar 
disks, where they increase the gas temperature 
through photoelectric heating and provide the sur- 
face for chemical reactions, a better characteri- 
zation of their chemical properties is needed to 
improve our understanding of circumstellar disk 
structure and evolution. Furthermore, the life cy- 
cle of hydrocarb on molecules in space is poorly 
understood (e.g.. iHenning fc Salamal T998I ) . Here, 
we aim to contribute to these topics. 

2. Data set 

2.1. Spitzer IRS spectra 

We have gathered infrared spectra of 53 iso- 
lated Herbig Ae star s obtained with the Infrared 
Spectrograph fIRS. ,Houck et al.ll2004l) aboard the 
Spitzer Space Telescope. Most of these sources 
were observed within the programs The miner- 
alogy of proto-planetary disks surrounding Her- 
big Ae/Be stars (PID 3470, PI J. Bouwman) and 
Probing the disk mineralogy and geometry of Her- 
big Ae/Be stars (PID 20308, PI B. Acke). The 
spectra are characterized by a very high signal- 



The 6.8-/xm feature is linked to CH bending modes 
in both CH2 and CH3 functional groups, the 7.2- 
fim feature only to CH3. Hydrocarbon molecules 
with mixed hybridization display aromatic as well 
as aliphatic emission features. 

The work presented in this paper is performed 
in the framework of an infrared study of a large 
sample of Herbig Ae stars, based on IR spec- 
tra collected wit h the Spitzer Space Telescope. 
Acke et al. ( 20091 ) investigated the geometry of the 



dust disk based on the strength and slope of the 
thermal infrared emission. We refine the interpre- 
tation of the Meeus group I/group II classification 
and show that the height of the inner disk varies 
over the sample. This leads to a varying amount 



^|20i^ noise ratio, typically several hundreds. For a 
detailed description of the data reduction process, 
the reduced spectra, and the sample targets, we 
refer to JIO. 

For each of the sample targets, published photo- 
metric data are gathered to comp ile the spectra l 
energy distributions. A reddened iKurucd (Il99lh 
model was fitted to the observed UV-optical pho- 
tometry. The circumstellar IR excess flux is char- 
acterized by several parameters such as excess 
magnitudes and spectral indices in the near-to- 
far-IR and at (sub-)mm wavelengths. Table [1] 
summarizes a few of these parameters. Based on 
the SED, we have classified the sample into the 
M eeus groups according to the method proposed 
by Ivan Boekel et al. ( 2005 ). For details, we refer 
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to JIO. 

2.2. Extraction of the PAH spectrum 

In this section, we describe how the hydrocar- 
bon emission features are separated from the ther- 
mal emission of the circumstehar dust. The IRS 
spectrum was split into two wavelength intervals 
(5-7 and 7-14 /im) which are treated differently. 

The 5-7-/im wavelength range is characterized 
by the presence of carbonaceous features on top of 
a relatively smooth dust continuum. In this inter- 
val, the continuum was approximated with a sin- 
gle spline through anchor points at 5.35, 5.45, 5.58, 
5.82, 6.66, 7.06, 7.37, 7.40, 7.55 and 7.70 /im. Note 
that the broad 7-9- fim PAH feature is treated as 
a continuum to the weaker features that are lo- 
cated on top of it. The features in the continuum- 
subtracted spectrum are measured according to 
the procedure described in Sect. 12.31 

The PAH features in the 7-14-/im wavelength 
range are less easily extracted from the original 
spectrum, as the underlying "continuum" may 
contain strong feat ures of amorphous and crys- 
talline silicate dust. iKeller et al.l (|2008l ) show that 
spline fits are only applicable to spectra with a 
high feature-to-continuum contrast. In spectra 
with stro ng silicate b a nds, artificial features are 
created. Keller et al. ( 20081) therefore only con- 
sider the strongest PAH sources in their analy- 
sis. In an attempt to reach weak PAH features, 
we follow a different approach. Six of our sample 
targets have a spectrum characterized by strong 
PAH emission and by the absence/weakness of 
silicate features. These sources are HD 34282, 
RR Tau, HD 97048, HD 135344B, HD 141569 and 
HD 169142. Their PAH profiles are representa- 
tive for those observed in the entire sample, cov- 
ering the observed range in feature peak posi- 
tion, shape and width. A single spline was fitted 
through continuum anchor points at 5.52, 5.83, 
6.66, 7.06, 9.54, 10.32, 13.08 and 13.85 /xm. The 
spline-subtracted spectra were used as PAH tem- 
plate spectra. Each template spectrum is cut into 
five separate "PAH feature complexes" at 5.5-7.0, 
7.0-8.2, 8.2-9.2, 10.5-12.3 and 12.2-13.5 /im. The 
procedure to obtain full 5-16-/im spectral decom- 
position of the IRS spectra makes use of these tem- 
plates. For each of the PAH complexes, the fit- 
ting routine picks one of the six template spectra, 
appropriately rescaled. The best fit for a certain 



target includes a PAH spectrum which is a scaled 
combination of five template spectra. The inclu- 
sion of PAH template spectra prevents the fitting 
routine from compensating for the PAH features 
using the dust species in the model. Hence, the re- 
sults of the dust fit should be more reliable. More- 
over, after subtraction of the dust-|-continuum 
part of the best spectral fit to the original spec- 
trum, the residual spectrum is the real PAH spec- 
trum of the target. Fig. [T]shows the residual spec- 
tra. A detailed description of the dust modeling 
can be found in JIO. 

The procedure works well for most spectra. 
However, some residuals of the crystalline sili- 
cate features remain. This is due to a slight mis- 
match in peak wavelength between the observa- 
tions and the computed silicate opacities (see JIO 
for details). The residuals hamper an accurate 
extraction of the PAH features with low peak- 
to-continuum ratios, especially in the 10-14-/im 
region. Fig. [2] shows the resulting PAH spectra 
of two silicate-dominated targets. In the case of 
HD 35187, the procedure leads to an excellent re- 
sult; for HD 244604, the extraction method cre- 
ates strong artifacts in the 7-12-/im region. The 
features in this region cannot be detected with 3- 
sigma significance. 

The aromatic-CCC bending-mode feature is 
present in the 16-19-/im range of a handful of 
sources. A spline through anchor points at 14.0, 
14.8, 15.0, 15.4, 15.6, 18.1, 18.6, 19.3 and 19.5 /im, 
approximating the underlying continuum, is sub- 
tracted. 

2.3. Derivation of the observables 

We aim at investigating whether the aliphatic 
and aromatic features observed in the spectra dis- 
play the signature of chemical diversity and evo- 
lution of their carriers. To do so, the detected 
features are captured in a number of representa- 
tive parameters. We stay as close as possible to 
the real measurements, and avoid the use of ap- 
proximate line profiles such as Gaussian curves. In 
this section, we describe how the observables are 
obtained. 

For each feature, a wavelength interval indi- 
cating its blue and red limits is defined based on 
the best-quality PAH-dominated spectra where 
the feature is detected (HD 34282, HD 97048, 
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HD 169142 and RR Tau). This wavelength inter- 
val was visually checked for each individual target 
in order to avoid false detections. To minimize the 
effect of residuals of the continuum subtraction, ei- 
ther the spline or dust model, a straight line which 
connects the blue and red limit is subtracted from 
the spectrum. The integrated line flux is the to- 
tal integral of the resulting continuum-subtracted 
feature. Furthermore, the centroid position of 
the feature (i.e. the flux- weighted barycenter), the 
full width at half maximum (FWHM), the peak 
position and the peak flux are determined. Uncer- 
tainties on the first two parameters were deduced 
directly from the error bars on the spectrum. For 
the FWHM, peak position and peak flux, the fol- 
lowing approach was taken: 100 spectra of the 
feature were simulated around the original spec- 
trum, assuming a Gaussian error distribution with 
a sigma equal to the error in each spectral pixel. 
In each of these spectra, we determine the FWHM, 
peak position and peak flux of the feature. The 
standard deviation of these measurements is taken 
as the error on the feature parameter determined 
from the original spectrum. 

To quantify the influence of our choice of the 
feature wavelength limits on the measurements, 
the procedure described above was repeated over 
a slightly larger wavelength interval: the blue limit 
was shifted one spectral pixel to the blue, and 
the red limit one pixel to the red. The difference 
between the two values for each parameter was 
used to estimate the error induced by the choice 
of wavelength limits. The final error on our mea- 
surements is the square root of the quadratic sum 
of the latter, systematic, error and the observa- 
tional error described above. 

Some of the features are a blend of individual 
features. When the spectral resolution allows it, 
the components of these complexes are separated. 
The largest feature is considered to be part of the 
underlying continuum of the smaller feature(s). In 
general, a spline is fitted to the wing of the large 
feature on one side of the smaller feature, and to 
the continuum (which is zero in the continuum- 
subtracted spectra) on the other side. The resid- 
ual spectrum is that of the smaller feature, the 
larger feature is set equal to the spline in the region 
where the small feature is located. The parame- 
ters of both features are determined as described 
above. Fig. |3] shows the procedure for the 6-/im 



region. This method was also used to separate the 
7-dll and 8.6-fim features and the 10.6, 11.3, 12.0 
and 12.7-fim features (Fig. HJ. Again, an a priori 
fixed set of anchor points was used for all stars in 
the sample. 

From the line fluxes of the detected aromatic 
features, i.e. all features except those at 6.8 and 
7.2 fim, a total PAH flux is computed. The same 
is done for the aliphatic features. The ratio of the 
aliphatic and aromatic line fluxes is equal to their 
luminosity ratio iaii/-^PAH- We also compute the 
feature-to-stellar luminosity ratio for each 

feature / separately. The stellar luminosity is that 
of the Kurucz model fitted to the SED. 

3. Analysis 

3.1. Detected features 

A detected feature is defined as a feature with a 
continuum-subtracted peak flux exceeding the 3- 
sigma level. Strong aromatic features are detected 
at the well-known wavelengths 6.2, 7.8, 8.6, 11.3 
and 12.7 /im. Secondary features are observed at 
5.7, 6.0, 10.6, 12.0 and 13.5 /xm. We did not detect 
the feature at 14.3 ^m, linked to the quintet CH 
out-of-plane bending mode, in any of the sample 
sources. The 16-19-//m complex is measured in 
five targets (see Fig. [5]) . Aliphatic features are de- 
tected in a number of sources, at 6.8 and 7.2 /im. 
In total, 13 spectrally separated features are iden- 
tified. 

In 37 of the 53 sample sources (70%), one or 
more significant PAH features have been detected. 
The detection rate is 16/20 (80%) in group I, and 
21/33 (64%) in group II. Aliphatic features are 
found in 22 of the 40 sources for which a short- 
wavelength spectrum is available (55%). The 
aliphatic detection rate is 11/15 in group I (73%) 
and 11/25 in group II (44%). 

Table[2]gives an overview of the sample-average 
parameters of the detected features. We indicate 
the possible carriers of each of the features as men- 
tioned in the literature. Tables |3] and |4] list the 
measured line fluxes and upper limits. Table [5] 
gives the fractional feature-to-stellar luminosities 
and PAH-to-aliphatic luminosities. 



Because almost all our sample targets belong to Class B', 
with the dominant sub-band at 7.8 fim, we hereafter refer 
to this feature as the 7.8-fim feature. 
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3.2. Correlations 

We searched for correlations between the fea- 
ture parameters, the stehar parameters and the 
disk parameters derived from the SED. To com- 
pare two parameters, the Kendall tau-rank corre- 
lation coefficient wa s computed, w hich is a non- 
parametric statistic (jKendalllll938l ). A correlation 
is marked statistically significant if its p-value is 
less than 1%, under the assumption that no cor- 
relation is present, before and after removal of the 
outliers. Hence, apparent correlations based on 
only a few data points are removed, as well as 
those induced by outliers. 

The sample targets are classified in the two 
Meeus groups. In all figures presented in this pa- 
per, group I sources are the black dots, group II 
the red dots. For all feature parameters, we 
checked whether their values are significantly dif- 
ferent in one group with respect to the other. We 
used the non-parametric Wilcoxon rank-sum test 
for assessing whether two independent samples 
of observations come from the same distribution 



( Mann k Whitnevlll947n 



4. Results 

4.1. The chemical diversity of PAHs in 
Herbig Ae stars 

The Peeters spectral classification is interpreted 
in terms of chemical diversity of the hydrocar- 
bon carriers of the features. The CC stretch- 
ing feature at 6.2 /im in astronomical spectra is 
shifted with respect to its theoretical position at 
6.4 /im. This is also the case in our sample. The 
shift has been ascribed to the inclusion of nitrogen 
atoms in the PAH skeletal structure. The feature 
is increasingly bluer when the N atom is substi- 
tuted deeper in the PAH molecule ( Hudgins et al.l 
20051). Also the anion-to-cati on ratio may play a 
role (jBauschlicher et allbOOOh . However, the peak 
positions of the 6.2 and features correlate 

(Fig. |6]). While the general blueshift of the fea- 
ture from 6.4 to 6.2 /zm may still be attributed to 
nitrogen pollution, the latter cannot account for 
the pr ogressive sh i ft of the 6.2 and 7.8-^m fea- 
tures. iPino et al.l (|2008l ) have performed labora- 
tory measurements and show that a shift of the 
6.2-fim feature to longer wavelengths points to a 
higher aliphatic/aromatic content of the hydrocar- 



bon mixture. 

The feature-to-stellar luminosity ratios of al- 
most ah PAH features (at 5.7, 6.0, 6.2, 7.8, 8.6, 
10.6, 11.3, 12.0 and 12.7/^m) correlate with each 
other. This indicates that these PAH features ap- 
pear together in the spectrum of a Herbig Ae star, 
and roughly scale together in strength. The rel- 
ative strengths of the aliphatic and aromatic fea- 
tures, on the other hand, are not connected. We 
also find dozens of significant correlations of line 
fiux ratios. The strongest correlation in our sam- 
ple is between the 6.2 and 7.8-/im features, both 
linked to CC modes. Tight connections were also 
found between the 6.2 and 8.6-/im, the 7.8 and 
8.6-/im, and the 11.3 and 12.7-/xm features. Fig. [7] 
shows two examples of such correlations. The 
6.2, 7.8 and 8.6-/xm features, attributed to ionized 
PAHs, are closely related to each other, as well as 
the features in the ll-14-/im range, due to out-of- 
plane CH bending modes of neutral and ionized 
PAHs. 

Honv et al.l (|200l[ ) and iKenereLaD (|2008h dis- 



cuss the link between the 6.2/11.3 and 12.7/11.3 
line fiux ratios. The first is interpreted as a mea- 
sure for the degree of ionization, the second as 
a measure for the irregularity of the PAH's edge. 
The authors conclude that both ratios point to the 
degree of processing of the molecules. An increase 
in PAH molecular size also enhances the 11.3-/im 
feature with respect to the 6.2 and 12.7-/xm fea- 
tures. We confirm this correlation between line 
flux ratios (p-value 1%). iKeller et all (|2008[ ) also 
flnd a decrease of the 12.7/11.3 ratio with decreas- 
ing stellar temperature. Our data do not show a 
clear correlation for the entire sample, but a corre- 
lation may be present for the group I sources only 
(Fig.|8]). In any case, it appears that the 12.7/11.3 
ratio is different in group I and II. We come back 
to this issue in the Sect. 14.41 



4.2. The 16-19-/im feature 

Five stars (HD 34282, HD 36917, HD 97048, 
HD 100453 and HD 135344B) display an emission 
feature at 16-19 fim. The complex is attributed to 
the out-of-plane skel etal modes of la rge elongated 
PAH molecules (Peeters et al.ll2004l) . The detec- 



tion is remarkable because our sample targets are 
Class B' sources, while this feature complex is 
mostly seen in ISM-like Class A objects. 
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The strength and shape of the 16-19-/zm fea- 
ture is not correlated to those of other features at 
shorter wavelengths. This is consistent with the 
interpretation that the short- and long- wavelength 
features are produced by different PAH popula- 
tions. Four of the five targets with detected 16- 
19-/im emission have a very red SED: their 30/13.5 
continuum fiux ratios are among the highest in the 
sample (see Table Possibly, these sources are 
still surrounded by an envelope or extended outer 
disk in which the chemistry of the PAH molecules 
is governed by the tenuous but harsh interstellar 
UV field rather than tha t of the central star (see 
also iBoersma et ah 2008 ). 

Alternatively, the detection or otherwise of the 
feature may be due to a contrast effect. The con- 
tinuum emission at 20 fim is dominated by warm 
dust grains close to the central star, while the fea- 
ture is produced in the outer disk via stochastic 
heating. The extreme 30/13.5 continuum flux ra- 
tio could indicate that, while the outer disk is still 
strongly flared, the inner disk has started to set- 
tle. An additional indication for this hypothesis 
could be the absence of the 10-/J,m silicate emis- 
sion feature in HD 34282, HD 97048, HD 100453 
and HD 135344E|1: when silicate dust coagulates 
and grows beyond several fims, the emission bands 
in the mid-IR range disappear. The settling of the 
inner disk implies a reduction of the warm con- 
tinuum emission at 20 fim and hence an increased 
feature-to-continuum contrast. Later on, when the 
outer disk eventually starts to settle as well, the 
strength of the 16~19-/im complex decreases to- 
gether with the degree of flaring. 

4.3. Hydrocarbon molecules and silicate 
grains 

Based on the IR spectra, the properties of the 
(sub-)/im-sized silicate grains in the upper layers 
of the disk can be derived. JIO show that only a 
few dust species contribute to the emission-band 
spectrum: amorphous silicates with olivine and 
pyroxene stoichiometry, the Mg-rich crystalline sil- 
icates forsterite and enstatite, and silica. Different 
grain sizes (0.1, 2 and 5 ^m) are included in the 



2 On the other hand, HD 37411, RR Tau, HD 141569 and 
HD 169142 are also PAH sources without a 10-/im siUcate 
feature, but their spectra do not show the 16-19-/^m fea- 
ture. 



spectral fit. The fitting routine was performed on 
two wavelength ranges, 5-17 /im and 17-35 /zm, to 
sample the hot dust close to the star and the cooler 
dust further out and deeper in the disk. This 
yielded the mass fraction of each of the species, 
for each grain size and in both intervals. For de- 
tails, we refer to JIO. 

We have taken the dust mass fractions and 
searched for correlations with the hydrocarbon 
feature parameters derived here. None were found 
with a p- value below 5%, with one marked excep- 
tion. The FWHM of the 6.2-/im feature appears to 
be strongly correlated to the mass fraction of the 
0.1-/im silica grains, derived from the spectral fit 
to the 17-35 /zm spectrum (Fig. [5]). Targets with a 
broad 6.2-/im feature have less small silica grains 
in the cold/outer disk. It is conceivable that this 
is a statistical false-positive correlation out of the 
~3000 investigated comparisons between silicate 
and PAH parameters. However, the correlation 
has a very low p- value (1 x 10^^). The probability 
that such a strong random correlation would occur 
here is only 3%. Also the fact that both parame- 
ters were measured in disjunct wavelength ranges 
in the spectrum, strengthens our confidence in the 
reality of the correlation. 

The relation is unexpected and we can only 
speculate about its origin. PAHs are stochasti- 
cally heated and their emission comes from fur- 
ther out in the disk surface, where grains in ther- 
mal equilibrium are much colder. It is likely that 
the PAHs producing the 6.2-/im feature and the 
cold silica grains are located in the same physi- 
cal region. (jHudgins et al.ll2005l ) mention that pe- 
ripheral oxygen groups can influence the 6.2-/im 
feature position. Perhaps the cumulative effect 
of a population of PAH molecules with such pe- 
ripheral groups broadens the 6.2-/im feature. The 
observed correlation would then indicate that the 
0.1-/im silica grains and the PAH molecules are 
chemically competing for oxygen atoms: the oxy- 
gen is either locked in peripheral groups of PAH 
molecules, or included in 0.1- fim silica grains. 

4.4. Influence of the disk geometry 

We checked for differences between Meeus 
group I and II sources. This classification is based 
on the thermal infrared emission seen in the SED, 
and thus indicates the geometry of the dust disk. 
Group I sources have a fiared dust disk, while the 
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dust disk in the group II members has flattened, 
probably as a resuh of dust grain growth and 
sedimentation to the midplane. 

In a number of papers (iMeeus et al. 



Acke fc van den Anckeij|2004 : Habart et al 



20011 
20041) 



it has been claimed that the strength of the PAH 
features depends on the disk geometry. A flared 
disk has a larger illuminated surface. Group I 
sources are therefore expec ted to produce stron ger 
PAH emission. However, iKeller et al.l (|2008l ) do 
not find this correlation in their sample of Herbig 
Ae/Be stars. Here, we elaborate on this. 



Ma king use of radiative transfer models. iMeiier et al 
()2008 ') demonstrate that the degree of flaring of 
the dust disk largely depends on one parameter: 
the mass in sub-/Ltm-sized dust grains. A higher 
dust mass results in a higher opacity at optical 
wave lengths, a warmer outer disk, and more flar- 
ing. Acke et al. (200^ show in their Fig. 3 that 
the flux excess above the stellar photosphere at 
60 /zm is a tracer for the dust mass, and hence for 
the degree of flaring. Fig. [TU] compares the IRAS 
60-/im excess to the excess at 30 /im. The latter 
is determined from the Spitzer IRS spectra and 
the Kurucz model for the staiH. A very strong 
correlation is noted, which indicates that also the 
excess at 30 /im traces the disk geometry. Because 
IRAS 60-/Ltm photometry is not available for all 
sample sources, and to rely only on the Spitzer 
IRS spectra, we use the excess at 30 fim in the 
following. 

We investigate whether the luminosity of indi- 
vidual PAH features, relative to the stellar lumi- 
nosity, correlates with the excess at 30 /zm. The 
former is the fraction of stellar flux captured and 
re-emitted by the PAH molecules. We find that 
the feature-to-stellar luminosity of most of the fea- 
tures (5.7, 6.0, 6.2, 7.8, 8.6, 10.6, and 12.7 ^m) 
increases significantly with increasing 30- //m ex- 
cess. Fig. [TT] shows the correlation for the most fre- 
quently detected feature, at 6.2 ^m. Also the total 
PAH-to-stellar luminosity ratio correlates with the 
far-IR excess (p- value 9 x 10~^). 

It is clear that the disk geometry influences 
the PAH-to-stellar luminosity ratio. However, the 
large amount of scatter indicates that the degree 
of flaring is not the only parameter. Source-to- 
source differences in PAH abundance are certainly 



important. Furthermore, there are targets which 
have a flattened dust disk (group II), but are still 
strong PAH emitters. A solution to this apparent 
discre pancy has been pro p osed in a f ew recent pa- 



;rep 

pers. Ivan der Plas et al.l (|2008l ) and iFedele et al 



(|2008l ) show that the dust disk can decouple from 
the gas disk, probably due to dust grain growth 
and subsequent settling to the midplane. This re- 
sults in a flattened dust disk, whi le the gas disk 



can st ill have a flared structure. iVerhoeff et al 
(|2010l) lave further elaborated on this idea and 
performed detailed modeling of HD 95881, a tar- 
get with an extremely blue group H-type SED, but 
strong PAH emission. 



Along the same lines, iDuUemond et al.l ([20071) 
have modeled dust sedimentation in protoplane- 
tary disks with PAHs. They flnd that, if the 
PAH abundance stays the same, sources with 
sedimented dust disks should display more pro- 
nounced PAH features than disks with perfect ver- 
tical mixing. Two effects are important: the re- 
duction of the continuum emission (increased con- 
trast) and the increased amount of UV-exposed 
PAH molecules (increased PAH luminosity). Un- 
der the assumptions of the models, a decrease in 
far-IR continuum emission should be accompanied 
by an increase of PAH lumino sity. The oppo- 



site tr end is observed, however. iDullemond et al 



^The excess magnitude is defined as M = 2.5 log _Ftot/^* 



(|20Q7l ) speculate that a signiflcant fraction of the 
PAHs is removed through aggregation. In any 
case, a reduction of PAH abundance in the flared 
gas disk above the dust disk, e.g. through PAH 
photodestruction, is necessary to reconcile theory 
and observations. 

We found indications that there is indeed a 
small difference in PAH chemistry between flared 
and flattened disks. There is a strong correlation 
between the peak flux ratios 11.3/6.2 and 11.3/7.8 
(Fig.[T2|). This is not surprising, as the 6.2 and 7.8- 
/im features are both CC modes and thus strongly 
linked to each other (see Sect. 14.11) . From the flg- 
ure, however, it seems that the distributions of 
the 11.3/6.2 and the 11.3/7.8 ratios are shifted to- 
wards lower values in group II. Also the 12.7/11.3 
line flux ratio appears to be higher in this group 
(0.29 ± 0.08 vs. 0.22 ± 0.05 in group I, p-value 
0.5%). This suggests that the 11.3-^m feature is 
weak in the targets with a small far-IR excess, and 
may indicate that the PAH molecules in targets 
with a flattened dust disk are more processed, i.e. 
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rough-edged and/or smaller, than those in flared 
disks. 

4.5. Aliphatic versus aromatic hydrocar- 
bons 

Almost all stars in our sample display Class B' 
PAH features, but it is clear that even within this 
class, spectral differences occur. Here we argue 
that these changes are mainly due to the stellar 
UV field to which the hydrocarbon molecules are 
exposed. 

The centroid (and peak) position of the 7.8-/im 
feature decreases with increasing effective temper- 
ature. Fig. [T3] shows the data for our sample of 
Herbig Ae stars. We have expanded the sample 
with a set of pre-main-sequence stars of different 
stellar masses, obtained in other Spitzer IRS pro- 
grams. On the low-mass end, we have included 
five T Tauri stars with detected PAH emission. 
RX J1842.9-3532, RX J1852.3-3700, HD 143006, 
RX J1612.6-1859Aand IRXS J132207.2-693812 
were observed within the Formation and Evolu- 
tion of Planetary Systems (FEPS) Legacy Sci- 
ence Program (PID 148, PI M. Meyer). The 
spectra an d stellar effective tem peratures were 
taken from iBouwman et al. ( 2008f) . On the high- 
mass end, six Herbig Be stars were added. The 
spectra of HD 290770, Hen 3-180, LkHa 257, 
PDS 216, PDS 344 and RNO 6 were acquired 
within the Spitzer Program Structure and compo- 
sition of disks surrounding Herbig Be stars (PID 
50180, PI M. van den Ancker). Their effective 
temperature has been estimated based on their 
spectral type. The correlation between 7.8-/im 
centroid position and Teg is clear for the Herbig 
Ae stars, and is even more prominent when in- 
cluding the lower- and higher-mass young stars. 
The contribution of the red sub-bands appears to 
be more significant in cooler stars than in hot- 
ter stars, which is supported by the fact that the 
width of the 7.8-/xm feature increases with redshift 
(Fig. [141). Finally, also the peak and centroid po- 
sition of the 6.2-/im feature are correlated to the 
stellar effective temperature (p- value 5 x 10~^ and 
4 X 10~^ resp.): cool stars display red features. 

The correlation between T^e and 7.8-/xm wave- 
length position was alread y found in smaller sam- 
ples by 'Sloan et al.' ('2007t). iBoersma etHI (|2008l ) 
and tKcller et al. (20081). As we have shown, it 
holds for a broad variety of pre-main-sequence 



stars, w ith spectral types ran ging from MO to 



r es ran ; 
(|2010l ) have recently ex- 
tended the correlation to lower Tes and redder 
features with the detection of hydrocarbon emis- 
sion in S-type AGB stars. The stellar radiation 
field strongly infiuences the chemistry of the sur- 
rounding hydrocarbon molecules, regardless of the 
object's evolutionary status or specificities of the 
circumstellar environment. 



Sloan et al. 



()2007() propose that a red feature 
indicates the presence of a hydrocarbon mixture 
with a higher aliphatic/aromatic ratio, an idea 
that is supporte d by laboratory measurements 



porte c 

(jPino et al.l l2008l ). Here we report a few newly 
found correlations in our sample of Herbig Ae stars 
that provide further evidence for this interpreta- 
tion. 

The position and width of the 7.8-/im feature 
are linked to the effective temperature of the cen- 
tral star. Moreover, the strength of the aromatic 
8.6-fim CH in-plane bending mode, relative to the 
6.2 and 7.8-/im CC features, and the aliphatic-to- 
aromatic luminosity ratic|f| depend on TcS as well. 
In fact, all these parameters are interconnected. 
Table IH] summarizes these correlations. Figs. [T51 
and [16] show that, while targets with a red 7.8- 
/im feature (i.e. cool stars) have a low fiux ratio of 
the aromatic 8.6-/im CH and 7.8-/im CC features, 
they produce fractionally more aliphatic emission. 
It is clear that illuminated hydrocarbon molecules 
around stars with different effective temperature 
are chemically different. 

Herbig Ae stars are by nature strong photo- 
spheric UV emitters. Because of the rather low 
accretion rates in most targets, the additional UV 
excess due to accretion energy dissipation is rela- 
tively unimportant for these stars. We can there- 
fore consider the effective temperature as a direct 
indicator of the stellar UV radiation field. 

Strong UV radiation appears to preferentially 
destroy or convert the aliphatic component in 
the disks. Photodestruction of the CH bonds 
is expected from labor atory experiments (e.g., 
Munoz Caro et al1l200l[ ). This effect is stronger 
for aliphatic molecules, as aromatic rings give 
more stability to the molecule by allowing distri- 



''This is the ratio of the summed line fluxes of the aliphatic 
6.8 and 7.2-fim features with respect to those of the de- 
tected PAH features in the covered wavelength range. 
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bution of the UV-induced excitation energy. This 
is consistent with the lower hne fluxes of the 8.6- 
lim feature relative to the ahphatic CH features 
observed in cooler stars. The disappearance of 
the 6.8 and 7.2-^m features, both linked to CH 
bonds in aliphatic hydrocarbons, could therefore 
point to a more rapid dehydrogenation of the 
aliphatic molecules, rather than a complete re- 
moval of this component. Howey er, the experi- 
ments of iMunoz Caro et al.l (|2001l ) show that the 
aromatic CH bonds are also eventually destroyed 
by the UV radiation. The aromatic feature at 
8.6 /im should therefore become weaker with in- 
creasing UV strength as well. The opposite trend 
is observed: sources with strong UV fields display 
a higher aromatic CH/CC flux ratio. We note that 
this interpretation relies on the identification of 
the 8.6-/im feature as due to aromatic CH modes. 
Alternatively, it may be that also aromatic CC 
modes contribute to its strength. In the latter 
case, the observations show that the aromatic hy- 
drocarbons undergo structural changes induced by 
UV field, resulting in Class A features and an in- 
creased contribution of CC modes to the 8.6-/im 
feature. 

The observed increase in the aromatic CH 
strength at higher te mperatures is consisten t with 
experimental work bv lMennella et ah ( 200l[ ). who 
suggest that the heating of carbon grains should 
reduce the number of aliphatic CH bonds and in- 
crease the number of aromatic CH bonds. Also 
UV irradiance increase the relative degree of aro- 
matic clustering, but in thermal annealing at tem- 
peratures above 300 K, this increase is accom- 
panied with a growth of the graphitic clusters. 
Under UV irradiation , the clusters remain small 



( Mennella et al.lll998[ ) 



In diffuse interstellar clouds, large (> 30 C 
atoms) PAHs are expected to have normal hy- 
drogen coverage (i.e. one H atom per peripheral 
C atom; aromatic bond), while very large PAHs 
may be fully hydrogenated (i.e. tw o H atoms; 
aliphatic bond, lie Page et alll2003l ). The tran- 
sition of full to normal hydrogen coverage could 
occur through UV photolysis of the aliphatic CH 
bonds and replacement with aromatic CH bonds. 
This would reduce the aliphatic fingerprint in the 
infrared spectrum and pump the aromatic CH/CC 
ratio. Because UV processing keeps the cluster 
sizes small, thermal annealing may be needed to 



grow the clusters and increase the aromatic CH 
emission strength. The observational trends in our 
sample then indicate that not only UV photode- 
struction, but also thermal annealing of carbon 
grains and hydrogenation in the UV-shielded disk 
interior, is important in circumstellar disks around 
Herbig Ae stars. 

5. Conclusions 

In this paper, we have investigated hydrocarbon 
emission in the Spitzer IRS spectra of a sample 
of 53 Herbig Ae stars. We find that 70% of the 
sources display PAH emission, and 55% aliphatic 
features at 6.8 and/or 7.2 /im. Five targets (9%) 
display a feature complex at 16-19 /im. Below we 
list our main conclusions. 

• We confirm the correlation between the PAH 
feature-to-stellar luminosity ratio and the 
shape of the infrared spectral energy dis- 
tribut ion, first mentioned by iMeeus et al 



POOII ). PAHs are a tracer of the gas disk, 
while the SED probes the thermal dust emis- 
sion. If the dust and gas disk are colocated, 
the trend indicates that the illuminated sur- 
face of a fiared disk is larger than that of 
a flattened disk. The volume occupied by 
excited PAHs is therefore larger in the first 
group. It is suggested in the literature that 
sources with settled dust disks can nonethe- 
less produce strong PAH emission if their 
gas disk is still flared and exposed to stellar 
UV flux. This could explain the existence 
of the few targets with flattened dust disks 
but strong PAH emission. However, most 
targets follow the general trend. We con- 
clude that PAHs in the dust-poor disk atmo- 
sphere are either rapidly removed from the 
gas phase, or that the gas and dust in most 
disks are well coupled and settle together. 

The stellar radiation field determines the 
chemistry of the hydrocarbon population in 
the circumstellar disk. Ranked according 
to increasing stellar effective temperature, 
the infrared spectra of Herbig Ae stars show 
an increase of the aromatic CH emission, 
but a decrease of the aliphatic CH emis- 
sion. Aliphatic CH bonds are more easily 
destroyed than aromatic CH bonds, which 
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explains the decrease in aliphatic signature. 
However, this effect cannot explain the in- 
crease in the aromatic CH/CC ratio with 
increasing UV strength. The solution may 
be found in thermal annealing and the hy- 
drogcnation balance of the PAHs. The lat- 
ter can only influence the PAH chemistry 
if turbulent vertical mixing is important at 
the disk's surface, and brings PAH molecules 
from UV-immersed regions to the disk inte- 
rior and back. 

• The 6.2 and 7.8-//m CC features are redder 
in cooler stars. The link with the relative 
strength of the aliphatic features indicates 
that the wavelength shift indeed is a measure 
for the aliphatic/ aromatic content ratio of 
the hydrocarbon mixture. 

• The strong connection between spectral fin- 
gerprint and stellar effective temperature 
shows that hydrocarbons immediately re- 
act to the UV field to which they are ex- 
posed. Therefore, the chemistry of hydrocar- 
bon molecules does not in general trace disk 
evolution. However, the observations indi- 
cate that targets with a fiattened dust disk 
harbor hydrocarbons which may be more 
processed. They appear to have rougher 
molecular edges and/or are smaller than 
those in fiared dust disks. 



the Jet Propulsion Laboratory, California Insti- 
tute of Technology under a contract with NASA. 
BA thanks Kees DuUemond for helpful discus- 
sions. 



We have expanded the knowledge of the chem- 
istry of hydrocarbon molecules in disks. Once 
irradiated, the free-flying molecules are quickly 
chemically altered by the stellar UV field and 
their history is deleted. However, PAHs in disks 
must have a significant aliphatic component when 
they are appear at the disk surface. Labora- 
tory experiments have shown that the formation 
of aliphatic (CH2) groups on the edges of PAH 
molecul es can already occ ur under very mild con- 
ditions (jjager et al.ll2006l T « 400 K). Such con- 
ditions are available in the warm regions of disks 
around Herbig Ae stars, which makes it probable 
that in-situ production and (re-)hydrogenation of 
the hydrocarbon mixture are processes which take 
place below the disk surface. 



This work is based on observations made with 
the Spitzer Space Telescope, which is operated by 
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Fig. 1. — Residual PAH spectra. Next to the source name, the Meeus group classification is indicated. As a 
reference, the PAH spectrum of HD 169142 is overplotted in red, scaled to the peak of the residual spectrum. 
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Fig. 2. — (a) Spitzer IRS spectrum of HD 35187 and (b) PAH spectrum, extracted by subtracting the best 
fit to the dust and continuum (JIO). The major features at 6.2, 7.8, 8.6, 11.3 and 12.7/im are clearly visible 
(dashed lines), as well as the aliphatic features at 6.8 and 7.2 /im (arrows), (c) Spitzer IRS spectrum of 
HD 244604 and (d) PAH spectrum. Residuals of the silicate dust fit are present on a few-percent level and 
hamper a reliable detection of the (weak) PAH features longward of 7 /im. 
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Fig. 3. — Extraction of the features in the 6-/im range of the IRS spectrum of HD 169142. A continuum 
sphne, through the anchor points indicated with +, was subtracted. The 5.7, 6.8 and 7.2-/xm features are 
not blended and can be measured directly. A sphne is fitted to the blue wing of the large 6.2-/j,m feature 
at well-chosen wavelength positions (diamonds) . It is assumed that the large feature reaches the continuum 
level (i.e. zero) at 6 /im. The resulting line profiles (colored lines) are plotted on top of the original spectrum 
(black). See text for details. 
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Fig. 4. — Extraction of the features in the 7-9-^m range (top) of the residual PAH spectrum of HD 169142. 
Sphnes are fitted at a priori fixed anchor wavelengths (diamonds). The resulting line profiles (colored lines) 
are plotted on top of the original spectrum (black). Note that the 7.2-/im feature is extracted, considering the 
7.8-/im feature as continuum emission (see also Fig. [3l). Bottom: Extraction of the features in the 10~13-/im 
region. 



19 



15.5 16.0 16.5 17.0 17.5 18.0 18.5 

Wavelength [fim] 



Fig. 5. — Detected features of the skeletal CCC out-of-plan e bending mode. Th e spectrum of the NGC 7023 



reflection nebula is shown at the bottom as a reference ( Sellgren et al.l 12007 ). It is rebinned to a lower 



resolution of 160 to allow direct comparison with the Herbig Ae spectra. The full vertical line indicates the 
position of the H2 0-0 S(l) molecular line. Even if it would be present, confusion with the broad PAH feature 
can be excluded. 
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Fig. 6. — The peak position of the 6.2 and 7.8-/im features (p-value 3 x 10 ^). The black dots represent 
group I sources, the red group II. This convention is used throughout the entire paper. 
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Fig. 7. — Strong correlation between the line flux ratios 6.2/6.0 and 7.8/6.0 (p-value 1 x 10 ^, top) and 
12.7/6.0 and 11.3/6.0 (p-value 2 x lO"'^, bottom). The 6.2 and 7.8 -/im features are both CC stretching 
modes, the 11.3 and 12.7-/ini features are due to CH out-of-plane bending modes. 
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Fig. 8. — Line flux ratio of tfie CH out-of-plane bending features at 12.7 (duo+trio CH) and 11.3-/im (solo 
CH) as a function of stellar effective temperature. A correlation including all sample targets can be excluded 
(p- value 15%), but is tentatively present for the group I sources alone (p- value 5%). 
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Fig. 9. — Correlation between the FWHM of the PAH 6.2-/im feature versus the mass fraction of the 0.1-/iin 
siUca grains in the cold dust component (p- value 1 x 10^^). 
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Fig. 10. — The excess magnitude at 30 and 60 /im. Both are strongly connected (p- value 2 x 10 ^). 
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Fig. 11. — The luminosity ratio of the PAH 6.2-/im CC stretch feature and the central star increases with 
increasing excess at 30-iim (p-value 4 x 10^^). 
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Fig. 12. — Correlation between the peak flux ratios 11.3/6.2 and 11.3/7.8 (p-value 2 x 10 ^). Note that the 
group II sources (red dots) have lower values than group I sources (black) for both ratios. 
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Fig. 13. — The centroid position of the 7.8-/xin feature decreases with in creasing stellar effective temperature. 



We have included the T Tauri stars with detected PAH emission from iBouwman et al.l (|2008l ). The Herbig 



Be stars were observed under Spitzer program PID 50180. The correlation is significant for the Herbig Ae 
stars alone (p- value 4 x lO^'^) and even more so for the entire sample of young stars (p- value 7 x 10^*). 
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Fig. 14. — A red 7.8-//m feature has a larger FWHM (p-value 2 x 10"^). 



29 



0.25 



00 

\. 

CD 
00 

o 



0.20 - 



0.15 



0.10 - 



0.05 



0.00 



-1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r- 



I II 1 



m 



i 



_l I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L_ 



7.8 



7.9 8.0 
Centroid position 7.8 [/^.m] 



8.1 



Fig. 15. — When the centroid wavelength of the 7.8-/im feature increases, the 8.6/7.8 hne flux ratio decreases 
(p- value 4 x lO""*). The 7.8-/zm feature is due to CC stretches, the 8.6-/ini feature is attributed to a CH 
in-plane bending mode in aromatic hydrocarbons. 
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Fig. 16. — Correlation between the aliphatic-to-PAH luminosity ratio and the centroid position of the 7.8-/im 
feature (p-value 4 x 10~^). 
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Table 1 
Stellar and disk parameters 





roup 


-L eff 


ivi ^_ 1 ^iii ) 




[13.5/7]'' 


[OUy J.O.OJ 








(mag) 


(mag) 








HD 31293 


I 


9520 


4.3±0.1 


9.3±0.2 


1.14±0, 


,05 


4.5±0.1 


HD 31648 


n 


8720 


4.6±0.1 


8.2±0.2 


1.28±0, 


,03 


1.19±0.03 


HD 34282 


I 


8720 


4.1±0.1 


9.3±0.2 


0.56±0, 


,09 


10.0±0.6 


HD 35187 


II 


8970 


2.8±0.1 


7.8±0.2 


2.5±0, 


,2 


2.1±0.1 


HD 35929 


II 


6870 


2.0±0.1 


3.9±0.2 


0.70±0, 


.03 


0.43±0.01 


HD 36112 


I 


7850 


3.6±0.1 


8.6±0.1 


1.37±0, 


.05 


4.1±0.2 


HD 244604 


II 


8730 


3.8±0.1 


7.9±0.2 


1.51±0, 


.05 


1.43±0.03 


HD 36917 


II 


10010 


2.8±0.1 


6.2±0.2 


0.86±0, 


.06 


1.35±0.06 


HD 37258 


II 


8970 


4.1±0.1 


7.9±0.2 


1.65±0, 


.04 


1.10±0.03 


BF Ori 


II 


8985 


4.2±0.1 


7.6±0.2 


1.36±0, 


.03 


0.87±0.03 


HD 37357 


II 


9230 


3.3±0.1 


7.8±0.1 


1.70±0, 


.06 


1.85±0.04 


HD 37411 


II 


9100 


4.1±0.1 


9.0±0.1 


1.12±0, 


.08 


4.0±0.2 


RR Tau 


II 


8460 


5.0±0.1 


8.9±0.1 


1.1±0, 


.1 


1.72±0.04 


HD 37806 


II 


9480 


5.1±0.1 


8.3±0.2 


1.25±0, 


.03 


0.85±0.07 


HD 38120 


I 


10500 


4.7±0.1 


10.5±0.2 


4.5±0, 


.5 


2.6±0.2 


HD 250550 


I 


10700 


5.0±0.1 


lO.liO.l 


2.2±0, 


.1 


2.5±0.1 


HD 259431 


I 


25400 




10.3±0.2 


- 




2.55±0.05 


V590 Mon 


I 


13000 


6.8±0.1 


12.2±0.1 


3.0±0, 


.6 


2.2±0.1 


HD 50138 


II 


12230 




14.4±0.2 


- 




0.78±0.04 


NX Pup 


II 


7290 




11.8±0.2 


- 




0.78±0.02 


HD 58647 


II 


10750 


3.6±0.1 


5.5±0.2 


0.66±0, 


.02 


0.49±0.01 


HD 72106S 


II 


9810 


3.4±0.1 


7.9±0.2 


2.0±0, 


.2 


1.62±0.07 


HD 85567 


II 


12450 


4.9±0.1 


7.9±0.2 


1.02±0, 


.02 


0.82±0.04 


HD 95881 


II 


8990 


4.9±0.1 


7.8±0.2 


l.OOibO, 


.03 


0.77±0.05 


HD 97048 


I 


10010 


3.9±0.1 


9.7±0.2 


1.6±0, 


,4 


5.9±0.4 


HD 98922 


II 


10500 




8.8±0.2 


- 




0.75±0.03 


HD 100453 


I 


7390 


3.9±0.1 


9.4±0.1 


1.7±0, 


1 


5.2±0.3 


HD 100546 


I 


10500 




11.3±0.2 


- 




3.5±0.2 


HD 101412 


II 


9960 


4.6±0.1 


7.9±0.2 


1.31±0, 


,04 


0.92±0.03 


HD 104237 


II 


8405 


3.4±0.1 


7.2±0.2 






1.28±0.03 


SS73 44 


I 


15000 


6.6±0.1 


11.5±0.1 


2.10±0, 


,09 


2.40±0.08 


HD 135344B 


I 


6590 


3.4±0.1 


7.9±0.2 


0.32±0, 


,02 


10.9±0.3 


HD 139614 


I 


7850 


3.0±0.1 


9.3±0.1 


4.1±0, 


,3 


4.2±0.3 


HD 141569 


II 


9520 


l.liO.l 


6.1±0.2 


0.73±0, 


,07 


6.8±0.2 


HD 142666 


II 


7580 


3.3±0.1 


7.4±0.2 


1.45±0, 


,06 


1.53±0.05 


HD 142527 


I 


6260 


3.1±0.1 


8.0±0.2 


0.95±0, 


,02 


5.0±0.1 


HD 144432 


II 


7345 


3.5±0.1 


8.0±0.2 


1.82±0, 


,06 


1.82±0.06 


HD 144668 


II 


7930 




8.3±0.2 






0.96±0.02 


HD 145263 


II 


7200 


0.8±0.2 


5.4±0.2 


1.9±0, 


,1 


2.0±0.1 


Wray 15-1484 


II 


30000 




9.0±0.2 






0.76±0.02 


HD 150193 


II 


8990 




9.1±0.2 






1.42±0.05 


HD 152404 


II 


6450 


2.5±0.1 


7.1±0.1 


1.13±0, 


,05 


3.3±0.1 


KK Oph 


II 


8030 




7.7±0.1 






1.04±0.02 


HD 158643 


II 


10100 


2.7±0.1 


5.2±0.2 


0.96±0, 


,02 


0.59±0.02 


HD 163296 


II 


8720 




6.4±0.2 






2.0±0.1 


HD 169142 


I 


8200 


2.8±0.1 


8.9±0.1 


1.7±0, 


,3 


7.8±0.5 


VV Ser 


II 


9000 


5.7±0.1 


8.6±0.2 


1.09±0. 


,02 


0.79±0.02 


T CrA 


I 


7200 




11.3±0.1 






5.0±0.3 


HD 179218 


I 


9810 




11.5±0.2 






2.4±0.2 


WW Vul 


I 


8430 


4.4±0.1 


8.8±0.1 


1.80±0, 


,05 


1.60±0.04 


HD 190073 


II 


8990 


4.3±0.1 


7.3±0.2 


1.14±0, 


,02 


0.75±0.02 


LkHa 224 


I 


7850 




8.9J-0,1 






1.99+0.02 


111) 20:!()_'l 


1 


8200 


:i,2 1 0.1 


!).! . 0,2 


5.(1 . 


() 


2.:i 1 0.2 



Note. — (a) The excess magnitude M(A) is relative to the stellar photosphere at A; (b) 
Continuum flux ratio at the given wavelengths. 
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Table 2 

Sample-averaged parameters of the detected features 



Feature 


Limits 


Mode 


Rcf. 


Det. 


<PP> 


<Cen> 


<FWHM> 


<logio 


L//L,> 


<log,„ L* /LpAH> 




(^tm) 






(%) 


(Aim) 


(^m) 


(^tm) 










5.7 


5.58-5.83 


CH-oop? 


1 


45 


5.69 ± 0.02 


5.70 ± 0.02 


0.07 ± 0.02 


-4.8 


± 0.6 


-2.2 ± 


0.3 


6.0 


5.83-6.10 


carbonyl-s? 


1 


78 


6.00 ± 0.04 


6.00 ± 0.03 


0.11 ± 0.02 


-4.2 


± 0.4 


-1.5 ± 


0.3 


6.2 


5.83-6.66 


CC-s (PANH?) 


1,5 


78 


6.26 ± 0.02 


6.28 ± 0.02 


0.16 ± 0.01 


-3.4 


± 0.5 


-0.7 ± 


0.2 


6.8 


6.66-7.06 


CHa and CH3 (ali) 


2 


38 


6.87 ± 0.03 


6.85 ± 0.07 


0.11 ± 0.04 


-4.7 


± 0.3 


-2.0 ± 


0.3 






or CH-ip n 


5 


















7.2 


7.06-7.36 


CH3 (ali) 


2 


50 


7.22 ± 0.02 


7.22 ± 0.01 


0.11 ± 0.03 


-4.6 


± 0.2 


-2.0 ± 


0.3 


7.8 


7.06-9.10 


CC-s 


1,5 


53 


7.88 ± 0.08 


7.94 ± 0.08 


0.7 ± 0.1 


-2.8 


± 0.6 


-0.24 ± 


0.06 


8.6 


8.25-9.10 


CH-ip 


1,5 


50 


8.66 ± 0.03 


8.64 ± 0.03 


0.25 ± 0.04 


-3.7 


± 0.6 


-1.1 ± 


0.2 


10.6 


10.35-10.90 


CH-oop s + 


3 


26 


10.64 ± 0.04 


10.61 ± 0.06 


0.25 ± 0.09 


-4.5 


± 0.7 


-2.1 ± 


0.3 


11.3 


10.35-12.30 


CH-oop s n( — ) 


3,5 


42 


11.27 ± 0.02 


11.3 ± 0.1 


0.29 ± 0.05 


-3.6 


± 0.5 


-1.0 ± 


0.2 


12.0 


11.70-12.30 


CH-oop d iiH — 


3 


28 


12.03 ± 0.05 


12.00 ± 0.03 


0.22 ± 0.08 


-4.5 


± 0.6 


-1.9 ± 


0.3 


12.7 


12.10-13.10 


CH-oop dt nH — 


3,5 


38 


12.7 ± 0.1 


12.61 ± 0.05 


0.41 ± 0.09 


-4.2 


± 0.6 


-1.7 ± 


0.1 


13.5 


13.20-13.90 


CH-oop q nH — 


3,5 


19 


13.59 ± 0.04 


13.6 ± 0.1 


0.19 ± 0.07 


-5.2 


± 0.5 


-2.7 ± 


0.4 


16-19 


15.60-18.10 


CCC-oop 


4 


9 


16.7 ± 0.4 


16.8 ± 0.3 


0.8 ± 0.3 


-4.0 


± 0.5 


-1.6 ± 


0.4 



Note. — The error bar represents the standard deviation in the sample. Det. — deteetion rate; PP — Peak position; Cen — Centroid position; FWHM 

— Full width at half maximum; hf/h^ — Feature-to-stellar luminosity ratio; Ly/LpAH — Feature-to-total- PAH luminosity ratio. The integration limits 
of the features are given in the seeond column. The third eolumn summarizes the identification of the mode of the feature, as published in the Uteraturc. 

OJ CH — carbon-hydrogen bond; CC — carbon-carbon bond, oop — out-of-plane bending mode; ip — in-plane bending or wagging mode; s — stretching 

mode. PANH — nitrogenated PAH; ali — aliphatic. For the features in the 10 — 14-/im region: s — solo or non-adjacent CH group; d — duo; t — trio; q 

— quartet CH group; n — neutral PAHs; -|- — cationic PAHs; — — anionic PAHS. 
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Table 3 

Line fluxes of the PAH and aliphatic features in the 5-10-/iM range 



Target 


5.7 


6.0 




6.2 


6.8 




7.2 


7.8 


8.6 


HD 31293 


<5 (- 


1.3) 


72±1 (- 


-13) 


475±4 (-13) 


24±4 (■ 


-1.3) 


15±2 (-13) 


181±6 (-12) 


24±6 (-12) 


HD 31648 


<7 (- 


13) 


224±9 ( 


-14) 


90±1 (-13) 


8±3 (- 


-13) 


92±8 (-14) 


<2 (-11) 


<5 (-12) 


HD 34282 


3±1 (- 


-13) 


96±3 (- 


-14) 


1010±5 (-14) 


17±4 (. 


-14) 


25±5 (-14) 


3380±8 (-14) 


387±6 (-14) 


HD 35187 


7±4 (- 


-14) 


88±4 (- 


-14) 


839±6 (-14) 


35±5 (■ 


-14) 


37±5 (-14) 


313±2 (-13) 


30±1 (-13) 


HD 35929 


<1 (- 


13) 


<7 (- 


14) 


<1 (-13) 


<4 (- 


13) 


<1 (-13) 


<5 (-13) 


<2 (-13) 


HD 36112 


12±9 (- 


-14) 


41±6 (- 


-14) 


414±7 (-14) 


36±7 ( 


-14) 


3±1 (-13) 


<1 (-11) 


<1 (-12) 


HD 244604 


<2 (- 


13) 


37±3 (- 


-14) 


148±3 (-14) 


<1 (- 


13) 


<8 (-14) 


<2 (-12) 


<9 (-13) 


HD 36917 


8±3 


-14) 


71±3 (- 


-14) 


364±4 (-14) 


14±5 ( 


-14) 


<4 (-13) 


123±1 (-13) 


20±1 (-13) 


HD 37258 


<2 (- 


13) 


<9 (- 


14) 


<1 (-13) 


<2 (- 


13) 


<9 (-14) 


<7 (-13) 


<5 (-13) 


BF Ori 


<9 (- 


14) 


llitl (- 


-14) 


61±3 (-14) 


<6 (- 


14) 


<8 (-14) 


<1 (-12) 


<4 (-13) 


HD 37357 


<9 (- 


14) 


24±2 (- 


-14) 


129±3 (-14) 


9±5 (- 


-14) 


<2 (-13) 


<3 (-12) 


<6 (-13) 


HD 37411 


11±2 (■ 


-14) 


24±3 (- 


-14) 


235±4 (-14) 


<9 (- 


14) 


<7 (-14) 


912±6 (-14) 


149±4 (-14) 


R.R, Tau 


48±7 (■ 


-14) 


94±3 (- 


-14) 


922±8 (-14) 


<1 (- 


13) 


<2 (-13) 


311±1 (-13) 


537±7 (-14) 


HD 37806 


<5 (- 


13) 


89±6 (- 


-14) 


479±9 (-14) 


<4 (- 


13) 


81±8 (-14) 


<1 (-11) 


<3 (-12) 


HD 38120 


<1 (- 


13) 


59±6 (- 


-14) 


412±8 (-14) 


19±7 ( 


-14) 


29±7 (-14) 


<7 (-12) 


<1 (-12) 


HD 250550 


9±3 (- 


-14) 


37±3 (- 


-14) 


166±3 (-14) 


15±8 ( 


-14) 


21±2 (-14) 


<6 (-12) 


<3 (-13) 


HD 259431 


n.s 




n.s. 




n.s. 


n.s 




n.s. 


n.s. 


n.s. 


V590 Mon 


<2 (- 


13) 


76±3 (- 


-14) 


681±4 (-14) 


<2 (- 


13) 


<3 (-13) 


243±2 (-13) 


31±1 (-13) 


HD 50138 


n.s. 


n.s. 




n.s. 


n.s. 


n.s. 


n.s. 


<8 (-12) 


NX Pup 


n.s. 


n.s. 




n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


HD 58647 


<9 (- 


13) 


<4 (- 


13) 


<5 (-13) 


<6 (- 


13) 


<6 (-13) 


83±2 (-13) 


12±1 (-13) 


HD 72106S 


29±5 (■ 


-14) 


144±3 ( 


-14) 


759±3 (-14) 


<2 (- 


13) 


19±6 (-14) 


162±3 (-13) 


30±2 (-13) 


HD 85567 


2±1 (- 


-13) 


76±6 (- 


-14) 


73±1 (-13) 


<5 (- 


13) 


43±7 (-14) 


232±4 (-13) 


34±3 (-13) 


HD 95881 


7±1 (- 


-13) 


163±6 ( 


-14) 


170±1 (-13) 


<6 (- 


13) 


6±2 (-13) 


569±5 (-13) 


83±3 (-13) 


HD 97048 


30±4 (■ 


-13) 


47±1 (- 


-13) 


551±3 (-13) 


15±3 (■ 


-13) 


7±5 (-13) 


2001±7 (-13) 


297±3 (-13) 


HD 98922 


n.s 




n.s. 




n.s. 


n.s 




n.s. 


n.s. 


21±1 (-12) 


HD 100453 


11±1 ( 


-13) 


39±1 (- 


-13) 


214±2 (-13) 


<9 (- 


•13) 


7±1 (-13) 


522±7 (-13) 


47±5 (-13) 


HD 100546 


n.s. 


n.s. 




n.s. 


n.s. 


n.s. 


n.s. 


23±2 (-12) 


HD 101412 


<2 (- 


13) 


56±3 (- 


-14) 


301±6 (-14) 


<5 (- 


13) 


2±1 (-13) 


137±2 (-13) 


27±1 (-13) 


HD 104237 


<6 (- 


13) 


<3 (- 


13) 


<4 (-13) 


<7 (- 


13) 


<4 (-13) 


<8 (-12) 


<6 (-12) 


SS73 44 


<1 (- 


13) 




13"! 


<^ f— 131 


<1 (- 


13) 


<9 ( — 14) 


f — 13^ 


<2 ( — 13) 


HD 135344B 


2±1 (- 


-13) 


11±1 (- 


-13) 


37±2 (-13) 


4±2 (- 


-13) 


33±7 (-14) 


1433±9 (-14) 


57±6 (-14) 


HD 139614 


<5 (- 


13) 


44±4 (- 


-14) 


712±5 (-14) 


4±2 (- 


-1.3) 


45±8 (-14) 


212±3 (-13) 


25±2 (-13) 


HD 141569 


11±5 (- 


-14) 


89±5 (- 


-14) 


524±6 (-14) 


<8 (- 


14) 


<1 (-13) 


1172±8 (-14) 


191±4 (-14) 


HD 142666 


20±8 (- 


-14) 


84±7 (- 


-14) 


73±1 (-13) 


73±7 (■ 


-14) 


43±6 (-14) 


245±3 (-13) 


25±2 (-13) 


HD 142527 


8±2 (- 


-13) 


339±8 ( 


-14) 


215±1 (-13) 


28±2 (■ 


-13) 


14±1 (-13) 


66±1 (-12) 


5±1 (-12) 


HD 144432 


16±7 (■ 


-14) 


47±5 (- 


-14) 


306±9 (-14) 


3±1 (- 


-13) 


42±5 (-14) 


<3 (-12) 


<2 (-12) 


HD 144668 


n.s. 


n.s. 




n.s. 


n.s. 


n.s. 


n.s. 


<2 (-12) 


HD 145263 


<2 (- 


13) 


<6 (- 


14) 


<1 (-13) 


<7(- 


14) 


<1 (-13) 


<1 (-12) 


<2 (-13) 


Wray 15-1484 


n.s. 


n.s. 




n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


HD 150193 


n.s. 


n.s. 




n.s. 


n.s. 


n.s. 


n.s. 


<9 (-12) 


HD 152404 


<2 (- 


13) 


15±4 (- 


-14) 


129±5 (-14) 


<2 (- 


13) 


<2 (-13) 


<1 (-12) 


<4 (-13) 


KK Oph 


n.s. 


n.s. 




n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


HD 158643 


<2 (- 


12) 


<i (- 


12) 


<1 (-12) 


<3 (- 


12) 


<8 (-13) 


<2 (-11) 


<3 (-12) 


HD 163296 


n.s. 


n.s. 




n.s. 


n.s 




n.s. 


n.s. 


<9 (-12) 


HD 169142 


10±2 (■ 


-13) 


442±9 ( 


-14) 


258±2 (-13) 


9±3 (- 


-13) 


67±9 (-14) 


729±2 (-13) 


89±1 (-13) 


VV Ser 


<5 (- 


13) 


25±4 (- 


-14) 


497±6 (-14) 


<1 (- 


13) 


<2 (-13) 


153±3 (-13) 


27±2 (-13) 


T CrA 


n.s. 




n.s. 




n.s. 


n.s 




n.s. 


n.s. 


<9 (-13) 


HD 179218 


n.s. 




n.s. 




n.s. 


n.s 




n.s. 


n.s. 


243±8 (-13) 


WW Vul 


<6 (- 


14) 


12±1 (- 


-14) 


79±1 (-14) 


<7(- 


44) 


<4 (-14) 


<6 (-13) 


<3 (-13) 


HD 190073 


<2 (- 


13) 


<2 (- 


13) 


<2 (-13) 


<5 (- 


13) 


<2 (-13) 


<3 (-12) 


<2 (-12) 


LkHa 224 


n.s. 


n.s. 




n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


HD 203024 


<2 (- 


13) 


<8 (- 


14) 


<1 (-13) 


<2 (- 


13) 


<1 (-13) 


<4 (-12) 


<1 (-12) 



Note. — The line fluxes and errors are in ergs cm ^ s ^ , with a (6) representing a X 10^. For non-detections, the 3-sigma upper 
hmit is given, n.s. = No spectrum available at this wavelength. 
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Table 4 

Line fluxes of the PAH features in the 10-14 ^um range. 



Target 


10.6 


11.3 


12.0 


12.7 


13.5 


HD 31293 


<5 ( 


-12) 


11±3 (-12) 


8±7 (- 


-13) 


<2 (-12) 


6±5 (- 


-13) 


HD 31648 


<2 ( 


-12) 


<7 (-12) 


<2 (- 


12) 


<2 (-12) 


<1 (- 


12) 


HD 34282 


39±2 


(-14) 


554±4 (-14) 


89±2 ( 


-14) 


138±2 (-14) 


5±2 (- 


-14) 


HD 35187 


11±8 


(-14) 


27±2 (-13) 


<3 (- 


13) 


103±6 (-14) 


<2 (- 


13) 


HD 35929 


<8 ( 


-14) 


<1 (-13) 


<7 (- 


14) 


<6 (-14) 


<6 (- 


14) 


HD 36112 


<7 ( 


-13) 


<2 (-12) 


<5 (- 


13) 


<4 (-13) 


<3 (- 


13) 


HD 244604 


<4 ( 


-13) 


<7 (-13) 


<4 (- 


13) 


<3 (-13) 


<3 (- 


13) 


HD 36917 


27±3 


(-14) 


534±6 (-14) 


59±3 ( 


-14) 


88±3 (-14) 


20±3 ( 


-14) 


HD 37258 


<2 ( 


-13) 


<9 (-13) 


<2 (- 


13) 


<2 (-13) 


<1 (- 


13) 


BF Ori 


<2 ( 


-13) 


<9 (-13) 


<2 (- 


13) 


<2 (-13) 


<1 (- 


13) 


HD 37357 


<3 ( 


-13) 


<2 (-12) 


<2 (- 


13) 


<2 (-13) 


<2 (- 


13) 


HD 37411 


<4 ( 


-14) 


65±7 (-14) 


<5 (- 


14) 


22±2 (-14) 


<5 (- 


14) 


RR Tau 


81±3 


(-14) 


658±7 (-14) 


41±2 ( 


-14) 


175±2 (-14) 


7±2 (- 


-14) 


HD 37806 


<1 ( 


-12) 


<9 (-12) 


<1 (- 


12) 


<1 (-12) 


<9 (- 


13) 


HD 38120 


<6 ( 


-13) 


<1 (-12) 


<5 (- 


13) 


<6 (-13) 


<4 (- 


13) 


HD 250550 


<1 ( 


-13) 


<9 (-13) 


<1 (- 


13) 


<3 (-13) 


<1 (- 


13) 


HD 259431 


34±1 


(-13) 


308±3 (-13) 


240 ±8 ( 


-14) 


73±1 (-13) 


44±7 ( 


-14) 


V590 Mon 


51±5 


(-14) 


29±1 (-13) 


<1 (- 


13) 


83±4 (-14) 


<1 (- 


13) 


HD 50138 


<3 ( 


-12) 


<3 (-11) 


<4 (- 


12) 


<1 (-11) 


<5 (- 


12) 


NX Pup 


<6 ( 


-13) 


<5 (-13) 


<2 (- 


13) 


<3 (-13) 


<2 (- 


13) 


HD 58647 


<2 ( 


-13) 


<3 (-13) 


<2 (- 


13) 


<1 (-13) 


<1 (- 


13) 


HD 72106S 


<3 ( 


-13) 


44±2 (-13) 


43±9 ( 


-14) 


120±8 (-14) 


<2 (- 


13) 


HD 85567 


<3 ( 


-13) 


39±2 (-13) 


5±1 (- 


-13) 


74±9 (-14) 


<2 (- 


13) 


HD 95881 


14±1 


(-13) 


90±2 (-13) 


8±1 (- 


-13) 


25±1 (-13) 


<3 (- 


13) 


HD 97048 


45±2 


(-13) 


633±4 (-13) 


57±1 ( 


-13) 


1547±9 (-14) 


18±1 ( 


-13) 


HD 98922 


59±5 


(-13) 


18±1 (-12) 


23±5 ( 


-13) 


75±8 (-13) 


<2 (- 


12) 


HD 100453 


<7 ( 


-13) 


95±4 (-13) 


<7 (- 


13) 


14±2 (-13) 


<5 (- 


13) 


HD 100546 


<5 ( 


-12) 


<6 (-11) 


<5 (- 


12) 


<1 (-11) 


<7 (- 


12) 


HD 101412 


<2 ( 


-13) 


151±7 (-14) 


<1 (- 


13) 


38±6 (-14) 


<1 (- 


13) 


HD 104237 


<3 ( 


-12) 


<9 (-12) 


<2 (- 


12) 


<2 (-12) 


<2 (- 


12) 


SS73 44 


<9 ( 


-14) 


<2 ( — 13) 


<8 (- 


14) 


<9 ( — 14) 


<9 (- 


14) 


HD 135344B 


<4 ( 


-14) 


260±3 (-14) 


48±1 ( 


-14) 


60±1 (-14) 


10±2 ( 


-14) 


HD 139614 


<3 ( 


-13) 


<5 (-13) 


<3 (- 


13) 


<1 (-12) 


<2 (- 


13) 


HD 141569 


12±2 


(-14) 


83±3 (-14) 


11±2 ( 


-14) 


29±1 (-14) 


4±1 (- 


-14) 


HD 142666 


2±1 


(-13) 


42±2 (-13) 


82±9 ( 


-14) 


11±1 (-13) 


<2 (- 


13) 


HD 142527 


<2 ( 


-12) 


161±8 (-13) 


<1 (- 


12) 


25±4 (-13) 


<1 (- 


12) 


HD 144432 


<1 { 


-12) 


<7 (-12) 


<9 (- 


13) 


<8 (-13) 


<7 (- 


13) 


HD 144668 


<7 ( 


-13) 


<9 (-12) 


<9 (- 


13) 


<1 (-12) 


<1 (- 


12) 


HD 145263 


<7 ( 


-14) 


<6 (-13) 


<6 (- 


14) 


<5 (-14) 


<5 (- 


14) 


Wray 15-1484 


33±4 


(-13) 


212±6 (-13) 


25±3 ( 


-13) 


52±3 (-13) 


6±2 (- 


-13) 


HD 150193 


<2 ( 


-12) 


<1 (-11) 


<2 (- 


12) 


<3 (-12) 


<3 (- 


12) 


HD 152404 


<3 ( 


-13) 


<2 (-12) 


<2 (- 


13) 


<1 (-13) 


<1 (- 


13) 


KK Oph 


<5 ( 


-13) 


<9 (-13) 


<4 (- 


13) 


<4 (-13) 


<3 (- 


13) 


HD 158643 


<1 ( 


-12) 


<2 (-12) 


<1 (- 


12) 


<1 (-12) 


<8 (- 


13) 


HD 163296 


<2 ( 


-12) 


<2 (-11) 


<2 (- 


12) 


<2 (-12) 


<2 (- 


12) 


HD 169142 


88±7 


(-14) 


208±1 (-13) 


206±6 ( 


-14) 


365±6 (-14) 


9±5 (- 


-14) 


VV Ser 


<3 ( 


-13) 


17±2 (-13) 


<3 (- 


13) 


<2 (-13) 


<2 (- 


13) 


T CrA 


<4 ( 


-13) 


<1 (-12) 


<5 (- 


13) 


<9 (-13) 


<7 (- 


13) 


HD 179218 


<4 ( 


-12) 


<3 (-11) 


<1 (- 


12) 


<7 (-12) 


<4 (- 


12) 


WW Vul 


<2 ( 


-13) 


<9 (-13) 


<1 (- 


13) 


<1 (-13) 


<9 (- 


14) 


HD 190073 


<9 ( 


-13) 


<2 (-12) 


<8 (- 


13) 


<7 (-13) 


<6 (- 


13) 


LkHa 224 


<6 ( 


-13) 


<7 (-13) 


<3 (- 


13) 


<5 (-13) 


<3 (- 


13) 


HD 203024 


<5 ( 


-13) 


<9 (-13) 


<4 (- 


13) 


<2 (-12) 


<3 (- 


13) 



Note. — ^The line fluxes and errors are in ergs cm ^ s ^, with a (6) representing a X lO**. For 
non-detections, the 3-sigma upper hmit is given. 
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Table 5 

pah-to-stellar and aliphatic-to-pah luminosity ratios 



Target 


Lpah/L* 


Lali/LpAH 


HD 31293 


44 ± 2 (- 


-4) 


13 ± 2 (-3) 


HD 31648 


25 ± 5 (- 


-4) 


28 ± 7 (-3) 


HD 34282 


1554 ± 5 ( 


-5) 


7± 1 (-3) 


HD 35187 


1228 ± 9 ( 


-6) 


15 ± 1 (-3) 


HD 36112 


13 ± 2 (- 


-4) 


28 ± 7 (-3) 


HD 244604 


13 ± 3 (- 


-4) 


<2 (-3) 


HD 36917 


602 ± 5 (- 


-6) 


5 ± 2 (-3) 


BF Ori 


171 ± 7 (- 


-6) 


<1 (-2) 


HD 37357 


9±1 (- 


4) 


10 ± 5 (-3) 


HD 37411 


34 ± 1 (- 


-4) 


<1 (-2) 


RR Tau 


2097 ± 8 ( 


-5) 


<2 (-2) 


HD 37806 


19 ± 4 (- 


-4) 


21 ± 5 (-3) 


HD 38120 


13 ± 3 (- 


-4) 


3 ± 1 (-2) 


HD 250550 


13 ± 1 (- 


-4) 


34 ± 9 (-3) 


V590 Mon 


328 ± 2 (- 


-4) 


<1 (-1) 


HD 58647 


83 ± 2 (- 


-6) 


<1 (-1) 


HD 72106S 


248 ± 3 (- 


-5) 


5 ± 2 (-3) 


HD 85567 


88 ± 1 (- 


-5) 


11 ± 2 (-3) 


HD 95881 


552 ± 4 (- 


-5) 


6±2 (-3) 


HD 97048 


901 ± 3 (- 


-5) 


6 ± 2 (-3) 


HD 100453 


47 ± 1 (- 


-4) 


8 ± 1 (-3) 


HD 101412 


190 ± 7 (- 


-5) 


11 ± 5 (-3) 


HD 135344B 


180 ± 2 (- 


-5) 


32 ± 8 (-3) 


HD 139614 


254 ± 4 (- 


-5) 


23 ± 6 (-3) 


HD 141569 


353 ± 2 (- 


-6) 


<1 (-1) 


HD 142666 


194 ± 2 (- 


-5) 


28 ± 2 (-3) 


HD 142527 


295 ± 5 (- 


-5) 


36 ± 2 (-3) 


HD 144432 


29 ± 1 (- 


-5) 


21 ±4 (-2) 


HD 152404 


106 ± 5 (- 


-6) 


<3 (-2) 


HD 169142 


605 ± 2 (- 


-5) 


11 ± 3 (-3) 


VV Ser 


411 ± 6 (- 


-5) 


<1 (-1) 


WW Vul 


347 ± 7 (- 


-6) 


<3 (-2) 



Note. — Only targets with a full 5-37-jum IRS 
spectrum and with detected PAH emission are 
included, a (6) represents a X lO''. 
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Table 6 
Correlations 



Parameters p-value 



Cen 7.8 /um 


Teff 




\ 


4 


X 


10- 


-3 


PP 7.8 /xm 


Teff 




\ 


2 


X 


10" 


-3 


FWHM 7.8 Atm 


Teff 




\ 


6 


X 


10- 


-3 


Cen 6.2 fj,m 


Teff 




\ 


4 


X 


10- 


-3 


PP 6.2 /xm 


Teff 




\ 


5 


X 


10- 


-4 


LF 8.6/6.2 


Teff 






2 


X 


10- 


-3 


LF 8.6/7.8 


Teff 




/ 


5 


X 


10- 


-3 


Lali/ LpAH 


Teff 




\ 


3 


X 


10- 


-2 


LF 8.6/ali 


Teff 






3 


X 


10- 


-3 


FWHM 7.8 iim 


Cen 7.^ 




/ 


1 


X 


10- 


-3 


LF 8.6/6.2 


Cen 7.1 




\ 


1 


X 


10- 


-3 


LF 8.6/7.8 


Cen 7.i 




\ 


4 


X 


10- 


-4 


Lali/ ipAH 


Cen 7.^ 




/ 


4 


X 


10- 


-3 


LF 8.6/ali 


Cen 7.i 


inm 


\ 


3 


X 


10- 


-3 


LF 8.6/6.2 


La.\i/ LpAH 


\ 


1 


X 


10- 


-2 


LF 8.6/7.8 


La\i/ LpAH 


\ 


4 


X 


10- 


-3 



Note. — A summary of the correlations indi- 
cating the dependence of the hydrocarbon chem- 
istry upon the stellar radiation field. The arrows 
indicate a positive or negative correlation. The 
p-valuc is the probability of a non-correlation. 
Cen: centroid position; PP: peak position; LF: 
line flux ratio; FWHM: full width at half maxi- 
mum; ali: the aliphatic 6.8 and 7.2-/im CH fea- 
tures. 
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